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Most studies of Coherent Synchrotron Radiation (CSR) have only considered the radiation from 
independent dipole magnets. However, in the damping rings of future linear colliders, a large fraction 
of the radiation power will be emitted in damping wigglers. In this paper, the longitudinal wakefield 
and impedance due to CSR in a wiggler are derived in the limit of a large wiggler parameter K. 
After an appropriate scaling, the results can be expressed in terms of universal functions, which are 
independent of K. Analytical asymptotic results are obtained for the wakefield in the limit of large 
and small distances, and for the impedance in the limit of small and high frequencies. 
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I. INTRODUCTION 

Many modern advanced accelerator projects Q, 0, El 
call for short bunches with low emittance and high peak 
current where coherent synchrotron radiation (CSR) ef- 
fects may play an important role. CSR is emitted at 
wavelengths longer than or comparable to the bunch 
length whenever the beam is deflected 0, . The strin- 
gent beam requirements needed for short wavelength 
SASE Free-Electron Lasers have led to intensive theo- 
retical and exp erimental studies IE IE El El EH EE 
EJ 0) El lid El over the past a few years where the 
focus has been on the magnetic bunch compressors re- 
quired to obtain the high peak currents. In addition 
to these single-pass cases, it is also possible that CSR 
might cause a microwave-like beam instability in stor- 
age rings. A theory of such an instability in a storage 
ring has been recently proposed in Ref. |l8| with exper- 
imental evidence published i n I l9l . Other experimental 
observations [2E 0, |5E HE l24j] may also be associated 
with a CSR-driven instability as supported by additional 
theoretical studies |2al2rl l27| . 

The previous study of the CSR induced instability 
assumed that the impedance is generated by the syn- 
chrotron radiation of the beam in the storage ring bend- 
ing magnets |Ts[ . In some cases (e.g. the NLC damping 
ring [28j), a ring will include magnetic wigglers which 
introduce an additional contribution to the radiation 
impedance. The analysis of the microwave instability 
in such a ring requires knowledge of the impedance of 
the synchrotron radiation in the wiggler. Although there 
have been earlier studies of the coherent radiation from a 
wiggler or undulator |23, y(| , the results of these papers 
cannot be used directly for the stability analysis. 
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In this paper, we derive the CSR wake and impedance 
for a wiggler. We focus our attention on the limit of 
a large wiggler parameter K because this is the most 
interesting case for practical applications. It also turns 
out that, in this limit, the results can be expressed in 
terms of universal functions of a single variable after an 
appropriate normalization. 

The paper is organized as follows. In Sec. 2, we write 
down equations for the energy loss of a beam in a wiggler. 
We then derive the synchrotron radiation wakefield in the 
limit of a large wiggler parameter K in Sec. 3. In Sec. 
4, we obtain the synchrotron radiation impedance for a 
wiggler, and in Sec. 5 we discuss our results. 



II. ENERGY LOSS AND LONGITUDINAL 
WAKE IN WIGGLER 

The longitudinal wake is directly related to the rate of 
energy loss d£ / dt of an electron in the beam propagating 
in a wiggler. For a planar wiggler, a general expression 
for d£ / dt as a function of the position s of the electron 
in the bunch and the coordinate z in the wiggler was 
derived in Ref. [3(j • We reproduce here the results of 
that work using the authors' notation: 

^-=e 2 k w f d S <D{s-s<,K,z)^±, (1) 
cdt J_ oa as' 

where A(s) is the bunch linear density, 

D(s,K,z) = \-2x 
s 

A - K 2 B(A, z) [sin A cos z + (1 - cos A) sin z] 

A 2 + K 2 B 2 (A,z) ^ ' 

B(A,z) = (1 - cos A - A sin A) cos z 

+ (A cos A — sin A) sin z , (3) 
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and A is the solution of the transcendental equation 

A / K 2 \ K 2 rr , A , A 

s = - ^l + _j +_{[2(l-cosA)-AsinA] 

x (cos A cos 2z + sin A sin 2z) - 2(1 - cos A)}. (4) 



In the above equations, we use the following dimension- 
less variables: s = 7 2 k w s and The parameter 
A is equal to k w (z—z r ), where z and z r are the projected 
coordinates on the wiggler axis of the current position of 
the test particle and the retarded position of the source 
particle, respectively. The internal coordinate s is defined 
so that the bunch head corresponds to a larger value of s 
than the tail. The wiggler parameter K is approximately 
K 93.4 B w X w , with B w the peak magnetic field of the 
wiggler in units of Tesla and X w the period in meters. In 
addition, 7 is the Lorentz factor, e is the electron charge, 
c is the speed of light in vacuum, and k w — 2ir/\ lu is 
the wiggler wavenumber. Note that the function D is 
a periodic function of z with a period equal to n. Also 
note that, desipte assuming K 3> 1, we still assume a 
small-angle orbit approximation, i.e., if/7 <C 1. 

We introduce the longitudinal wake W(s) of the bunch 
as the rate of the energy change averaged over the z co- 
ordinate: 



1 d£ f s 

w ( s ) = — 2-1: = ~ k ™ / ds'G{s 
e z cdt J_ x 



ds' 



where 



G(s) 



dzD(s,K, 



(6) 



w(s) = — fc. 



dG{s) 
ds 



so that 



W(s) 



ds'w(s-s')\(s'). 



(J) 



(8) 



x (cos Acos2f + sin A sin 2z) 
- 2(1 -cos A)}. 



In this limit, the expression for D, Eq. © can also be 
simplified: 

„. „ sin Acosz + (1 — cos A) s'mz , , 

m,z) = 2 B(A, z) (10) 

as long as A is not too small, A ^> 1/K. Again, the 
parameter K is eliminated from this equation. A detailed 
analysis supporting this approximation can be found in 
Appendix ^ 



III. WAKEFIELD 

Using Eq. © and JTJJ we find 

2 f* ..sin A cos z + (1 — cos A) sin I 

&(0--y o dz , (ii) 

where A = A(£, z) is implicitly determined by Eq. 
(JSJ. The integrand in this equation has singularities at 
points z = z s where Z?(A(£, z 8 ), z s ) — 0. It is shown 
in Appendix [5] that in the vicinity of a singular point 
B(A(C z)) cx (z— z^ 1 / 3 , and the singularity is integrable. 

We plot in Fig. ^the function G(Q calculated by nu- 
merical integration. A characteristic feature of the func- 
tion G is the presence of cusp points, at which the func- 
tion reaches local maxima and minima. An approximate 



and we dropped K from the list of arguments of the func- 
tion G. The positive values of W correspond to the en- 
ergy loss and the negative values imply the energy gain. 
The usual longitudinal wake w(s) corresponding to the 
interaction of two particles is then defined as 



Note that the wake Eq. (jJJ) is localized in front of the 
particle and vanishes behind it, w = for s < 0. 

In the limit of large K, we can neglect unity in the 
first bracket of Eq. |@J, assuming that K 2 /2 ^> 1. Such 
an approximation is valid, if we are not interested in the 
very short distances of order of (K k w (0.5A for 

the NLC damping ring wiggler |28|). We also introduce 
a new variable £ = s/ K 2 which eliminates the parameter 
K from Eq. |gj: 

C(A,z) - ^ + ^{[2(1 -cos A) -A sin A] 



-0.2 



-0.4 



CD 




(9) 



FIG. 1: The solid curve represents the G(() defined in Eq. 
lllll as a function of the normalized coordinate 2£/tt. The 
(x)-signs are the approximation given in Eq. 11 2> . 

location of these cusp points and the value of the func- 
tion G at these points can be understood with a simple 
physical argument presented in Appendix[n] It turns out 
that the minima are located at distances s between the 
particles equal to the integer number of the fundamental 
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radiation wavelength in the wiggler, and the maxima ap- 
proximately correspond to the distance equal to an odd 



number of half-wavelength. A simple analytical calcula- 
tion in Appendix O gives the following results 



J 



G(C) = 



4(2k+1)7t 



for C, = with n = 1,2, • 

fnr r ~ ( 2 "+ 1 ) 7r _ 1 
iU1 ~ 4 (2n+l)7r 



with n = 0, 1, • 



(12) 



These are the "x" points in Fig. ^ showing very good 
agreement with the numerical result. 



In Eq. we keep only the largest term 



A. Short-range limit 



B(A,z) = -A sin(A - z) 



(16) 



In the limit £ < 1, it follows from Eq. that A < 
1 as well. Eq. © can then be solved using a Taylor 
expansion of the right-hand side: 



24 C 



1/3 



(13) 



Expanding the integrand in Eq. Ijll|). keeping only the 
first non-vanishing term in A and substituting A from 
Eq. (JTJ) yields 



G{Q 



1 



^ (3C) 1/3 Jo 

43 2 /3r(±i) 



cos 2 / 3 , 



-1/3 



-0.99 C 



-1/3 



(14) 



The above result can also be obtained if one considers 
a wiggler as a sequence of bending magnets with the 
bending radius R = 7 /k w K\ cosi|. Indeed, in a bending 
dipole, the corresponding Gbcnd(s) = -2 s _1 / 3 /(3 i? 2 ) 1 / 3 
[JOI- Averaging Gbcnd over the wiggler period yields 
Eq. 114(1 . The reason why such a model gives the correct 
result in this limit, is that the formation length of the ra- 
diation is much shorter than the wiggler period, and one 
can use a local approximation of the bending magnet for 
the wake. 



B. Long-range limit 

In the limit £ ^> 1, the parameter A is also large, and 
Eq. can be further simplified: 



A sin A cos(A - 2z) 



(15) 



For D, one now finds, 



D(C,z) = 



with 



2 sin(z- A(C,z)) 2 ' 



(17) 



(18) 



where the function A(£, z) is implicitly determined by 
Eq. ((15JI . Averaging over one wiggler period, we find 



(19) 



with 



HO 



1 r 

- / dzF((, 



2tt 



smz 



dz . — 

sm(z — A) 



(20) 



It is easy to check that the function F is periodic, F((+ 
tt/2) = F(C), and F(0) = 0, F(tt/4) = -1 in agreement 
with Eq. fT5f> . The average value F(Q is equal to —1/2. 
Since F is periodic in £ with a period of 7r/2, using Eq. 
(I20II . we get a Fourier series representation for F(Q: 



no 



2n+ 1 



J, 



n+l 



2n + l 



cos(4(2n + 1)C) , 



(21) 
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where J n (x) is the Bessel function of the first kind. 
Derivations of the Fourier coefficients are presented in 
Appendix^! In Fig. El we plot the function F(C) defined 

I 



1 1 °° 



It is worth noting, that the asymptotic expression in 
the limit £ 1 in Ref. [3(J is incorrect — instead of the F- 
function the authors obtained a sine function, which only 
corresponds to the fundamental mode of the radiation 
and neglects contribution from higher-order harmonics. 



0.0 




FIG. 2: Plot of F(C) of Eq. (2J , 



The longitudinal wake defined in Eq. is plotted in 
Fig. 

IV. IMPEDANCE 

The impedance Z(k) is defined as the Fourier trans- 
form of the wake, 

/>oo 

Z{k) = / dsw{s)e~ lks 
Jo 

= -ik\ dCG(C)e- 4 ^ C , (23) 
7 Jo 

where ko = A^/ 2 k w /K 2 is the wiggler fundamental radia- 
tion wavenumber. 

We evaluated the integral in Eq. Q23|) using numeri- 
cally calculated values of the function G(£) in the inter- 
val [Cmin,Cmax], where Cmin ~ 10" 3 and (max ~ 50. The 



in Eq. I|21|) for one period. 

The corresponding long-range wake is then 



2 

cos(4(2n + 1)C) ■ (22) 
I 



8 
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FIG. 3: The normalized wake Green function 
—w(s)[K/('yk w )] 2 as a function of the normalized coor- 
dinate 2C/7T. 

contribution to the integral outside of this interval was 
calculated using asymptotic representations Eqs. I|14|) 
and (25) . 

The resulting imaginary and real parts of the 
impedance are shown in Fig. 01 and Fig. [S] respectively. 

The real part of the impedance can be related to the 
wiggler radiation spectrum I{u>) |3l| : 

ReZ(w) = ^ I(u) . (24) 

The spectrum I(lu) in the limit if > 1 is calculated in 
Appendix[E] It shows a perfect agreement with the result 
presented in Fig. [5] 

Simple analytical formulae for the impedance can be 
obtained in the limit of low and high frequencies. 

The low-frequency impedance corresponds to the first 
term in Eq. (55) for function G which does not oscillate 
with £: 

G(Q = ~. (25) 



Jn 



2n + 1 



Jn 



2n+ 1 



5 
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FIG. 4: The imaginary part of the normalized impedance 
Z(k)/k w as a function of the normalized wavenumber k/ko. 
Solid line — numerical solution from Eq. I23H . dotted line — 
analytical low-frequency asymptotic behavior from Eq. 1261 . 
and dashed line — analytical high-frequency asymptotic be- 
havior from Eq. 1271 . 
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Normalized wavenumber 

FIG. 5: The real part of the normalized impedance Z(k)/k w 
as a function of the normalized wavenumber k/ko ■ Solid line — 
numerical solution from Eq. 1231 , dotted line — analytical low- 
frequency asymptotic behavior from Eq. 1261 . and dashed 
line — analytical high-frequency asymptotic behavior from Eq. 



Using the definition in Eq. 12: il) , we then obtain the low- 
frequency asymptotic behavior of the impedance as 



Z(k) 



k 
«'o 



7e + log 



7T fc, 



k 
' k 



\ 2i 1 




1 log 


Go) 


7T 





(26) 



where, 7e ~ 0.5772 is the Euler Gamma constant. This 
asymptotic low- frequency impedance is plotted in Figs. 
0] and [S] for comparison with the numerical solution. 

Since we have an analytical expression for the short- 
range G(C) in Eq. I|14|) . we get the asymptotic high- 
frequency impedance as 



Z(k) 



6£[g 



A 



0.7HA 



(Kk u 



V 7 



7 

2/3 



2/3 



fc 1 / 3 



fc 1 / 3 



(27) 



with A = 3~ 1/3 r(2/3)(V3i-l) ~ 1.63i-0.94 18]. This 
asymptotic high-frequency impedance is plotted in Figs. 
^ and [S] for comparison with the numerical solution. 



V. DISCUSSION AND CONCLUSION 

In this paper, we derived the wakefield and the 
impedance for wigglers with K 2 /2 3> 1 due to the syn- 
chrotron radiation. Analytical asymptotic results are 



obtained for the wakes in the limit of small and large 
distances, and for the impedance in the limit of small 
and high frequencies. The results obtained in this pa- 
per are used for the beam instability study due to the 
synchrotron radiation in wigglers |32( . 
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APPENDIX A: DETAILS FOR DERIVING EQ. 

m 

In the limit of A < 1, according to Eq. Q, B — A 2 . 
In the numerator of the second term on the right hand 
side of Eq. we would have K 2 B(A, £)[sin A cos z + 
(1 — cos A) sin z ] - K 2 A 3 > A, as long as A » 1/K. 
This is allowed, since we are interested in the limit of 
K ^> 1. Hence, A is neglected in the numerator. In the 
denominator of the second term, then K 2 B 2 ~ K 2 A A ^> 
A 2 , as long as A > 1/K, hence A 2 is dropped. There- 
fore, the second term of Eq. @ is on the order of 1/A. 
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According to Eq. Q, in the limit of A < 1, we have 
s <~ K 2 A, so the first term of Eq. (J2J is on the order of 
1/(K 2 A), hence is much smaller than the second term in 
the limit of K ^> 1. Therefore the first term 1/s could 
be dropped. All these considerations lead us to Eq. (|10|l . 

For A ~ 1 and K ^> 1, according to Eq. we have 
s ~ X 2 > 1. Eq. © suggests that B(A,z) ~ 1. Now, 
in Eq. in the limit of K 3> 1, we can neglect A 

in comparison with K 2 in the numerator and A 2 in the 
denominator of the second term on the right hand side. 
We then note that the second term is on the order of 1, 
and is much larger than the first term 1/s ~ K~ 2 <C 1, 
hence we can drop 1/s to obtain Eq. QlOfl. 

Now let us study the limit of A ^> 1. Eq. J2J suggests 
that B(A,z) ~ A. For K > 1, then in Eq. ©, A and 
A 2 could be dropped in the numerator and denominator 
of the second term on the right hand side, respectively. 
The second term is on the order of 1/A. Now, according 
to Eq. |@J, in the limit of A > 1, we have a ~ K 2 A. 
Hence, in the limit of K ^> 1, the first term of Eq. J2J, 
which is on the order of l/(K 2 A), is negligible, compared 
with the second term. Hence, we obtained Eq. I|1U|) . 

So, in general, for large K, as long as A is not too 
small, i.e., A ^> 1/K, the simplification leading to Eq. 
ifTnf is always acceptable. 



APPENDIX B: SINGULAR POINTS IN D((,z) 

To find the scaling of the singularity, we assume that 
at the vicinity of the zeroes z s of B(A,z), the leading 
term scales as 



then we have 

B' ^ab(z- z s ) a ~\ (B2) 

where the prime indicates the derivative with respect to 
z. 

Let us first calculate B' . From Eq. Q we have, 



B' = — sin z + sin(z — A) 

+ cos(z - A)A - cos(z - A)AA' . 

To find A', we revert to Eq. ©, where, we find 

C(A,z) 



A' 



B 2 (A,z) 



(B3) 



(B4) 



with 



OA. : i = 2 sin (2.: - A) sin [ — ) A 



2 sin I j J - cos ) A 



(B5) 



Note that C(A, z) is a well defined function at the zeroes 
of B(A,z). 

Combining Eqs. (JbTJ, l|B3|) . (|B4l) . (|B5|) . we have, near 
the zeroes z s , 



B'(A,z) = 



b 2 (z-z s ) 2a ' 



(B6) 



B^b(z-z s ) a , 



(Bl) where, 



3 = - cos(z - A)A x C(A, ; 



\z=z B ;A=A(C,z B ) 



- cos(z - A)A |2 sin(2z - A) sin 



A x 



2 sin 



A 



— cos 



(B7) 



=i s ;A=A(C,z a ) 



Here, z s is defined as the solution of i?(A(£, z s ), z s ) = 0. 
Therefore, combining Eqs. (|B6(1 and i|B2|) . we find the 
scaling index a = 1/3 and b = (3 S) 1 / 3 . This means that 
D(C, z) has only an integrable singularity at z — z s with 
D oc (z - z s ) _1/3 . 



As a numeric illustration of the origin of the singular- 
ity, we plot in Fig. lU functions B(A((, z), z) and A(£, z) 
for C = 1.0. This plot shows that at the point where 
B = 0, both derivatives B' and A' are infinite, in accor- 
dance with Eqs. I|B2|) and JB4I) . 



APPENDIX C: SIMPLE PHYSICS MODEL 

Here, we give some explanation about the peaks and 
the zeroes in Fig. phased on a simple physics model. 

The CSR wake is actually the field emitted by a trail- 
ing particle which acts on the particle in front. Using a 
model presented in Ref. @ , the longitudinal force on the 
leading particle can be thought of as the component of 
the trailing particle's transverse Coulomb field projected 
onto the leading particle's direction of motion: 

W\\ = eEx(z) sin# , (CI) 

where Ej_(z) is the magnitude of the transverse Coulomb 
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_6 l ■ ■ ■ ■ ■ - 

12 3 

Z 



FIG. 6: A concrete example for ( = 1.0. Solid line: S(A(C = 
1.0, z), z) as a function of z; Dashed line: A(£ = 1.0, z) as a 
function of z. 




electric field at the retarded position from the trailing 
particle at retarded time. The argument z indicates that 
the amplitude of the transverse electric field is actually 
varying along the trajectory. In Eq. (|C1(I . is the angle 
between the direction of motion of the trailing particle 
and that of the particle at front. To illustrate the model, 
we give a schematic plot in Fig. \7\ 

To understand the wiggler wakefield, let us look at the 
four electrons in Fig. [S] The pair with solid arrows is 
separated by integer number times of the wiggler funda- 
mental radiation wavelength. During their journey, when 
the light emitted by the trailing electron cathches the 
electron in front, the instantaneous direction of motion 
of the front electron is always parallel to the direction 
of motion of the trailing electron at the retarded time 
when it emitted the light. Hence we have 9 = 0. So 
according to Eq. (|C1I) . the longitudinal force is always 
zero. This explains the zeroes in the longitudinal wake 
potential plotted in Fig. QJ. The pair with dashed ar- 
rows is separated by odd integer number times of half of 
the wiggler fundamental radiation wavelength. Averaged 
over one period, they make the largest angle between the 
instantaneous direction of motion of the front electron 
and the direction of motion of the trailing electron at the 
retarded time. Hence according to Eq. i|Cl(l . the longi- 
tudinal force reaches maximum. This explains the peaks 
shown in Fig. QJ. 

Let us calculate the values at these cusps points. Ac- 
cording to Eq. ©, when A = 2mr, we have £ = nir/2. 
According to Eq. Qllfl. the numerator of the integrand 
is zero at A = 2nir, while -B(A, z) = 2mtsmz accord- 
ing to Eq. (J3J). Hence the integral is zero. So we have 
G(nw/2) = for n = 1, 2, • • • . We find that it is true in 
the numerical solution in Fig. ^ 

At A = (2n + l)7r, according to Eqs. ijTT)) and ©, we 

get 



FIG. 7: Physics model to explain the longitudinal CSR wake. 




G(C) = - 



dz 



2 sinz 



2 cos z — (2n + l)w sin z 



4(2n+ 1)tt 
'4+[(2n+l)7r] 



for A = (2n + 1) 7T. (C2) 



According to Eq. ©, we have 



A 2 . A 1 

C = -; r cos z w - 

4 A 4 A 



for A = (2n+l)7r. (C3) 



FIG. 8: Physics model to explain the peaks and zeroes of 
G(s) shown in Fig. The solid line stands for the electron 
trajectory. The ellipses stand for the electrons. The arrows 
stand for the instantaneous direction of the motion. We group 
the four electrons into two pairs, one with solid arrows and 
the other with dashed arrows. 



In the above approximation, we average C, over one period 
in z. This becomes a good approximation, when A is 
large. This manifests itself in Fig. As we find from 
Fig. ^ with the increasing of £, therefore, the increasing 
of A, the above approximate value gets closer and closer 
to the numerical solution. Combining the results at the 
zeroes and the peaks, we get Eq. l JT2j l. 
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APPENDIX D: CALCULATION FOR THE 
FOURIER COEFFICIENTS 

Since we find that the function F(() is periodic in ( 
with a period of it/2, we could represent it in a Fourier 



series. The calculation for the Fourier coefficients are 
straightforward. We here illustrate one example. For 
to = 1,2, • • • , 



(F(C)cos[4(2m+ 1)C]) = 



d(F(() cos[4(2m + 1)C] 



\ f 2 dCcos[4(2TO+l)C] dz— Sm Z 



sin(z — A) 



•7T 2 , 

t Jo 



dA 



dz sin z sin(z — A) cos 



(2to+ 1)A 



(2to+ l)sin[2(A 



(2to+ l)sin2z" 



(Dl) 



Notice that, we have used the definition of F in Eq. Q2U[I. 
We also changed integral variable pair ((,2) to (A, z), 
using Jacobian obtained from Eq. I|15(l . To complete 
the integral in Eq. IjDIJI , we make use of the well-known 
identities 



cos(zcos#) = ^ e 2n (-) n J2n(z) cos2n8 , (D2) 

n=0 

with, eo — 1; £ n = 2 for n = 1, 2, • • • , and 

oo 

sin (z cos 

n=0 

All the other Fourier coefficients, including the average 
value (F(()) = —1/2, are obtained in the same manner. 
We then obtain the Fourier series representation in Eq. 



APPENDIX E: REAL PART OF THE 
IMPEDANCE AND THE WIGGLER RADIATION 
SPECTRUM 

Our approach to calculate I(u>) is based on the paper 
of Alferov et al. and that of Krinsky et al. |34j . 

The only difference is that we are dealing with large K 
case, hence we could eliminate K from the equations as 
was done for the impedance. Let us illustrate this in the 
following. The energy radiated per electron per unit solid 
angle per unit frequency interval per unit length is given 
by 



dI{Uj) -^±G m (K,,eA)H m (^) , (El) 

m— 1 x 



dCl 



where to is the harmonic number of frequencies in the 
spectrum. N is the number of periods of the wiggler. 
The polar coordinates 9 and cj) are defined so that, 8 = 



corresponding to the forward direction along the wiggler 
axis, and <j> = to be the plane of electron motion. For a 
large K , ui » 2 c k w r y 2 /(K 2 /2 + 7 2 9 2 ) is the fundamen- 
tal radiation frequency, and 



4 TO 



(i/Y 2 + 7 2 e 2 y 



{[Si 76* cost/) 



2 \ I K 2 + 7 2 « 2 

TO J 2j0 COS 
2 o2 „;„2 



+ ( 7 ^) 2 ^f sin" 



where, we define 



and 



and in turn 



if 2 



4(^+ 7 2 2 ) 



2i^ 7 6» cos<^ 

^- + 1 2 e 2 



n-\-m 



n— — oo 



and 



#2 = nJnim^) J 2n + m (m^ x ) 

n— — oo 

For finite, but large iV, 



sin 2 



AT 



7T 2 iV 2 ( ^ - TO 



(E2) 



(E3) 



(E4) 



(E5) 



(E6) 



(E7) 



determines the bandwidth of the radiation. In our theory, 
we focus on an infinite long wiggler, i.e., we are interested 
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in the limit of N — > oo. Under such condition, we have 
H m — ► u>i S(uj — mu)\)/N . Now we are ready to get the 
radiation spectrum by integrating over the entire solid 
angle. 



dO sin6> 



(E8) 



Due to the (^-function in H m , the integral over 9 could 
be done first. Also due to the fact that 9 is small, we 
approximate sin (9 m 9. Then the integral in Eq. (|E8|I is 
reduced into a 1-D integral as the following, 



A 2 00 

'<*>-&£ 

m— 1 



d(j) [ m — — 
k 



x < 



-I 2 



5*i cos < 
Sf sin 2 4>) 



Si 



2 cos < 



Notice, the integrand has a period of 7T, hence we need 
only integrate for one period. Si and $2 are further sim- 
plified as 



00 

E • 

n= — oo 



Si = 

X J2n+ 



k 




and 



5 2 = 



x J. 




2n+m 



(E10) 



(Ell) 



Now based on Eqs. (|2~4"|) and (|E9|) . we could compute 
(E9) the real part of the impedance as 



Re[Z(*)] = 2k n 



m—l 



dcj) [ m — — I x 
k , 



Si cos ( 



Si H S2 

m 1 2 cos< 



S 12 sin 2 ( 



(E12) 



Notice if is eliminated from this equation, as long as 
K is large. We sum up the first 30 harmonics to obtain 
the real part of the normalized impedance up to the 4th 
harmonic. Adding higher harmonics will not change the 



impedance within this frequency region within numerical 
accuracy. The result is identical to Fig. [SI supporting 
the validity of our calculation. 
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